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Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal 
freestream Mach number of 6 and Reynolds number of Re T « 460 are conducted at two 
wall temperatures ( T w /T r = 0.25, 0.76) to investigate the generated pressure fluctuations 
and their dependence on wall temperature. Simulations indicate that the influence of 
wall temperature on pressure fluctuations is largely limited to the near-wall region, with 
the characteristics of wall-pressure fluctuations showing a strong temperature dependence. 
Wall temperature has little influence on the propagation speed of the freestream pressure 
signal. The freestream radiation intensity compares well between wall-temperature cases 
when normalized by the local wall shear; the propagation speed of the freestream pressure 
signal and the orientation of the radiation wave front show little dependence on the wall 
temperature. 


Nomenclature 

C p heat capacity at constant pressure, J/(K-kg) 

C pp Space-time correlation coefficient of the pressure field, dimensionless 
C v heat capacity at constant volume, J/(K-kg) 

H shape factor, H = 5* /0, dimensionless 
M Mach number, dimensionless 

M r relative Mach number, M r = ( Uoo — U s )/a 00 , dimensionless 
Pr Prandtl number, Pr = 0.71, dimensionless 
R ideal gas constant, R = 287, J/(K-kg) 

Reg Reynolds number based on momentum thickness and freestream viscosity, Reg = <> ”^°° 9 , dimensionless 
Res 2 Reynolds number based on momentum thickness and wall viscosity, Res 2 = Poc / f°° e , dimensionless 
Re r Reynolds number based on shear velocity and wall viscosity, Re T = pwUrS , dimensionless 
T temperature, K 

T r recovery temperature, T r = T ao ( 1 + 0.9 * M , £,), K 

Ub bulk convection speed of pressure- carrying eddies, m/s 
U s convection speed of ‘effective’ radiation sources, m/s 
Uoo freestream velocity, m/s 

a speed of sound, m/s 

p pressure, Pa 

q dynamic pressure, Pa 

s entropy, J/K 

u streamwise velocity, m/s 
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u T friction velocity, m/s 

v spanwise velocity, m/s 

w wall- normal velocity, m/s 

x streamwise direction of the right-hand Cartesian coordinate 
y spanwise direction of the right-hand Cartesian coordinate 
z wall-normal direction of the right-hand Cartesian coordinate 
z T viscous length, z T = u w /u T , nr 
7 specific heat ratio, 7 = C p /C v , dimensionless 
S boundary layer thickness, nr 

S* displacement thickness, m 
k thermal conductivity, k = pC p /Pr , W/(m-K) 

9 momentum thickness, m 

p, dynamic viscosity, p = 1.458 x 10 ~~ 6 T ^ 0 4 , kg/(nr-s) 
v kinematic viscosity, v = p/ p, nr 2 -s 

p density, kg/m 3 

u> frequency, rad/s 

O vorticity, s _1 


Subscripts 

i inflow station for the domain of direct numerical simulations 

rms root mean square 

w wall variables 

00 freestream variables 

t stagnation quantities 


Superscripts 

+ inner wall units 

(•) averaged variables 

(•) / perturbation from averaged variable 

I. Introduction 

Prediction of laminar-turbulent transition is a critical part of the design of hypersonic vehicles because of 
the large increase in skin-friction drag and surface heating associated with the onset of transition. Despite 
continued advances in transition research, the physics of boundary layer transition over these vehicles is 
not fully understood due to the lack of detailed experimental measurements. Transition testing in conven- 
tional (i.e., noisy) wind tunnels has been an important avenue to understanding the transition behavior of 
hypersonic vehicles, despite the common knowledge that conventional wind-tunnel facilities cannot reliably 
simulate the in-flight transition behavior over a smooth surface due to the effects of the elevated levels 
of freestream disturbances . 1-3 Recent evidence suggests that freestream disturbances may also influence 
the accelerated onset of transition caused by isolated roughness elements on a nominally smooth surface . 4 
Transition measurements in low disturbance (i.e., quiet) wind tunnels better mimic the in-flight transition 
characteristics. However, because of the size and Reynolds number limitations of the existing quiet facili- 
ties, conventional tunnels will continue to be employed for the testing and evaluation of hypersonic vehicles, 
especially during ground tests involving large-scale models. 

Facility disturbances in conventional tunnels can impact not only the transition location but, possibly, the 
transition mechanism as well. As a result, the existing methodology to extrapolate wind-tunnel transition 
results to flight is rather crude and requires substantial improvement. To enable more effective use of 
the transition data from conventional facilities and permit more accurate extrapolation of the wind-tunnel 
results to flight, an in-depth knowledge of the broadband disturbance environment in those facilities must 
be developed. 

In unheated tunnels with adequate flow conditioning, the acoustic disturbances are likely to dominate 
the overall disturbance environment at Mach numbers of 2.5 or above , 5 ” 7 and their effect on transition 
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cannot be quantified in terms of a single metric corresponding to the root-mean-square amplitude of the 
freestream disturbances as indicated by the measurements at Purdue University. 3 * With the exception of 
the early measurements of freestream pressure fluctuations by Laufer 5 and a few others, 2,8 there are few 
measurements that are detailed enough to be suitable for either comparing with computational predictions 
or for developing models that can be used towards more reliable transition models. The measurements 
are typically susceptible to experimental errors due to the poor spatial resolution and/or limited frequency 
response of pressure transducers. Theoretical models for acoustic radiation from a supersonic boundary 
layer were developed by Phillips 10 and Ffowcs- Williams and Maidanik, 11 which attributed a major cause of 
the acoustic radiation to eddy Mach waves from boundary-layer turbulence converting supersonically with 
respect to the freestream. However, a lack of adequate knowledge concerning the boundary-layer turbulence 
restricted the theoretical predictions to the intensity of the freestream acoustic fluctuations alone. 

Direct Numerical Simulation (DNS) is a valuable tool that can overcome some of the aforementioned 
difficulties with both experimental measurements and theory and, hence, provide access to both flow and 
acoustic quantities that are difficult to obtain otherwise. The DNS can also isolate the acoustic radiation 
due to individual physical mechanisms, thereby avoiding any contamination due to secondary sources such 
as vortical and entropy fluctuations in the incoming stream. 

The successful application of DNS in capturing the freestream acoustic pressure fluctuations has been 
demonstrated for Mach 2.5 and Mach 6 boundary layers by the present authors. 12-14 Single and multi-point 
statistics of the computed freestream pressure fluctuations were reported, including the intensity, frequency 
spectra, space-time correlations, and convection velocities. However, a comprehensive understanding of 
the freestream disturbance field and its dependence of boundary-layer parameters (e.g., freestream Mach 
number, wall temperature, and Reynolds number) is still lacking. The current study focuses on exploring 
the dependence of boundary-layer induced pressure fluctuations on wall temperature at hypersonic speeds. A 
new DNS dataset at Mach 6 with a different wall temperature than the previous Mach 6 data 13 is introduced 
for the study of wall-temperature effects. 

The paper is organized as follows. The flow conditions selected for the numerical simulation of hypersonic 
turbulent boundary layers and the numerical method used for these simulations are outlined in Section II. 
Section III is focused on an analysis of freestream pressure fluctuations, highlighting their dependence on 
the wall temperature, conclusions thus far are outlined in Section IV. 

II. Simulation details 

A. Flow conditions and numerical methodology 

Table 1 outlines the freestream flow condition for the present simulations. The freestream condition is similar 
to the operational conditions of the Boeing/AFOSR Mach-6 Quiet Tunnel 3,4 (BAM6QT) under noisy flow 
conditions with a total pressure of 921 kPa and a total temperature of 433 K. Table 2 lists the values of 
the mean boundary layer parameters at the selected downstream location for statistical analysis ( x a = 5 7 <5, 
for Case M6Tw076 and x a = 80<5; for Case M6Tw025). The wall temperature for Case M6Tw076 is similar 
to that at the nozzle wall of BAM6QT, corresponding to a wall temperature ratio of T w /T r ss 0.76. Case 
M6Tw025 has the same freestream conditions and Reynolds number, Re T , as Case M6Tw025 but a lower 
wall temperature ( T w /T r ss 0.25). Thus, by comparing the results of Cases M6Tw025 and M6Tw076, the 
effect of wall cooling on the pressure fluctuations can be investigated at a fixed Reynolds number Re T . When 
necessary, the results of the two Mach 6 cases will be compared with a Mach 2.5 DNS case 14 to highlight 
the effect of freestream Mach number. 

The details of the DNS methodology, including numerical methods, initial and boundary conditions, 
have been documented in our previous papers. 13,14 Figure 1 shows the general computational set-up for 
Case M6Tw025. The computational set-up of Case M6Tw076 parallels that of Case M6Tw025. The inflow 
boundary condition is prescribed by means of a recycling-rescaling method and the selected rescaling length 
is sufficiently large to accommodate the eddy decorrelation length and minimize any spurious correlation in 
boundary layer fluctuations as a result of the recycling process. 14 Figures 3a and 3b show that the pressure 
fluctuations both at the wall and in the freestream for Case M6Tw025 have become homogeneous in the 
streamwise direction after x/Si = 55, indicating that both the boundary-layer turbulence and the freestream 
acoustic field have gone through the transient adjustment due to recycled inflow and reestablished equilibrium 
at the selected downstream location for statistical analysis. The streamwise variations of relevant quantities 
such as Cf, S , and rms pressure fluctuations at the wall and in the free stream are also used to gauge the 
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extent of the transient region. 

In the following section, averages are first calculated over a streamwise window of [x a — 0.9 5i,x a + 
0.9<5,:] (x a = 805.; for Case M6Tw025 and x a = 55(5; for Case M6Tw076) and spanwise locations for each 
instantaneous flow field; then, an ensemble average over flow field snapshots spanning a time interval of 
approximately 240(5;/I/ oo is calculated. To monitor the statistical convergence, flow statistics are computed 
by averaging over the whole or half the number of the flow-field snapshots and negligible difference (< 1%) 
is observed between the two. 

Power spectra are calculated using the Welch method 15 with 8 segments and 50% overlap. A Hamming 
window is used for weighting the data prior to the fast Fourier transform (FFT) processing. The sampling 
frequency is approximately 631/00/(5;, or 4 MHz, and the length of an individual segment is approximately 
53.2(5;/C/oo- Negligible difference in spectral estimation is observed within the reported frequency range, 
when the overall time record is subdivided in 12 segments instead of the baseline number of 8 segments. In 
addition, the DNS methodology has been extensively validated in previous work for supersonic/hypersonic 
turbulent boundary layers 16-19 and for supersonic shock wave-turbulent boundary layer interactions. 20-22 

B. Domain/grid sensitivity assessment 

In this section, the dependence of numerical results on the domain size and grid resolution for Case M6Tw076 
is investigated by a series of auxiliary DNS simulations listed in Table 3 where Case M6Tw076 is the baseline 
simulation and cases AI and All are two auxiliary DNS simulations for Case M6Tw076. All cases are identical 
in terms of numerical method, computational set-up, and flow conditions except for the parameters listed 
in the table. In particular, Case AI has the same grid resolution as the baseline case but only two-fifths of 
the spanwise domain size and 15% longer streamwise domain size, and Case All has the same domain size 
as case AI but with 1.5 times higher grid resolution in each direction when compared with the other two 
cases. As a result, the effects of streamwise and spanwise domain sizes can be investigated by comparing the 
results of the baseline case and AI, and the grid sensitivity can be demonstrated by comparing the differences 
between cases AI and AIL Furthermore, the effectiveness of the outflow boundary condition can be assessed 
by comparing the solutions near the outflow boundary for the baseline case and case AI. 

Figure 4a plots the intensity of pressure fluctuation across the boundary layer for DNS solutions cor- 
responding to each of these three cases. All curves collapse to within 1%, indicating that the single point 
statistics of most interest are insensitive to the spanwise domain size and the grid resolution. To further 
ascertain the dependence of high-frequency spectral components on the grid resolution, Figures 4b and 4c 
compare the frequency spectrum of pressure fluctuations at the wall and in the free stream, respectively, for 
the three simulations at the station selected for statistical analysis (Table 2). Good agreement is achieved 
up to ujS/Uoo ~ 25 or equivalently for u>v w /v%. « 1. The under-resolved frequency range (loS/Uoo > 25) 
contains less than 0.4% of the total energy and the pressure spectra in compensated form (not shown here) 
for the three cases are indistinguishable from each other over the range of resolved frequencies. 

To further assess the accuracy of the numerical scheme for resolving the propagation of acoustic distur- 
bances in the free stream, we note that the maximum frequency resolved by the streamwise grid spacing of 
Ax (note that the Ax > Ay > A z) may also be estimated as: 

_2nU c 

Umax ~ nAx [ ’ 

where U c refers to the streamwise propagation speed of the acoustic disturbance, and n denotes the minimum 
number of points per wavelength required for an algorithm to provide desired accuracy. Based on previous 
findings that the employed Weighted Essentially Non-Oscillatory (WENO) scheme can achieve a resolution of 
n = 8 PPW, 23 and setting UJU^ = 0.7 based on the results presented in Section D, the spatial resolution is 
deemed to be adequate up to u) max S /Uoo « 25. The latter bound approximately coincides with the frequency 
at which the differences between cases AI and All begin to appear. 

Similar assessments for domain size and the grid resolution are conducted for Case M6Tw025, and good 
frequency bandwidth is achieved up to oj max S /Uoo ~ 25. 14 Additional evidence to demonstrate satisfactory 
numerical accuracy of the current simulation is presented in the following sections, wherein the computed flow 
statistics are compared with both experiments and other high-quality simulations at similar flow conditions. 
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C. Comparison with experiments 

In this section, we demonstrate the accuracy of the current simulation by comparing the computed flow 
statistics with experiments at similar flow conditions. Figure 5 shows the comparison of DNS results (Case 
M6Tw076) with the wind-tunnel measurement and the calculation using Harris Boundary-layer code 24 con- 
ducted by Katya Casper for a Mach 5.8 turbulent boundary layer on the nozzle wall of BAM6QT under 
noisy-flow conditions (Re = 9.69 x 10 6 /m, P t>0 0 = 965 kPa, T t>00 = 429 K). The DNS and experiments agree 
well with each other in terms of both boundary-layer profile and wall-pressure spectrum. In comparison, the 
Mach number profile based on the boundary layer code exhibits larger differences from the measurement and 
the DNS in the outer part of the boundary layer. Moreover, Figure 5c and Figure 5d show that DNS suc- 
cessfully extends the measured spectra to higher frequencies. The resolution of the high-frequency region as 
well as the acoustic radiation due to these high-frequency fluctuations are especially important for studying 
the receptivity process associated with second-mode waves in hypersonic wind tunnels. 

Additional comparisons of DNS results with both experiments and other high-quality simulations are 
presented in the following sections. 


III. Results 

In this section, the effects of wall temperature on the pressure fluctuations induced by hypersonic turbu- 
lent boundary layers are investigated by comparing results between Case M6Tw025 and Case M6Tw076. 

A. Boundary layer statistics 

Next, turbulence statistics across the boundary layer are plotted and compared with well-established empir- 
ical correlations and numerical data sets. 

The van Driest transformed mean velocity profile based on the DNS is shown in Figure 6a. The mean 
velocity conforms well to the incompressible law-of-the-wall upon van Driest transformation and shows a 
(narrow) logarithmic region. In addition, the mean-velocity profile compares well with the experimental 
results by Bookey et al. 25 at Mach 2.9, Re T = 501 and the recent DNS results by Pirozzoli and Bernardini 26 
at Mach 2, Re T = 497. 

Figure 7 shows that wall cooling has a significant influence on p' rms /T w in the near-wall region ( z + < 150), 
with p' w rrns /T w « 3.4 for Case M6Tw025 and p' w , rrns / T w ~ 2.6 for Case M6Tw076, indicating an increase 
of approximately 30 percent in the cold-wall case. The pressure fluctuations have a rapid gradient near the 
wall for the cold-wall case, so that a portion of the pressure fluctuations damps out across the buffer layer. 
The influence of wall cooling becomes much weaker in the outer part of the boundary layer (z + > 150) and 
nearly vanishes in the free stream. Outside the boundary layer, p' rms /T w approaches a constant value of 
about 0.86 for the both Mach 6 cases. 

In terms of the influence of freestream Mach number, p' rms /T w shows only a weak Mach-number de- 
pendence within the boundary layer, consistent with previous findings. 13,2 ' However, it exhibits a strong 
Mach-number dependence in the free stream. The variation of freestream value with Mach number is con- 
sistent with the trend predicted by the experimental data reported by Laufer 5 (Figure 8). 

Similar to the pressure rms, the rms values of fluctuating velocity components (■ u' rms , v' rms , v' rms ) and 
other thermodynamic quantities (p' rms , T' rms ) also reach non-zero, nearly constant values within the free 
stream (at least up to z = 5(5 before the computational grid becomes significantly coarser). The constant 
asymptotes within the freestream indicate that local homogeneity in the wall-normal direction is established 
for all freestream disturbances, at least in terms of their rms amplitudes. 

B. Characteristics of freestream fluctuations 

The characteristics of freestream fluctuations will now be analyzed using the theory of modal analysis, which 
was initially proposed by Kovasznay 28 and further developed by Logan 29 and Smits and Dussauge. 30 Ac- 
cording to Kovasznay, 28 the fluctuations at any point within a uniform mean flow can be represented as a 
superposition of three different modes with co-varying physical properties: the vorticity mode, the acoustic 
or sound-wave mode and the entropy mode (also referred to as entropy spottiness or temperature spottiness) . 
Modal analysis has been used by a number of experimental researchers to characterize the freestream dis- 
turbance environment in their facilities. Examples include Donaldson and Coulter 2 who evaluated the level 
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and the sources of disturbances in the Arnold Engineering Development Center (AEDC) Tunnels A and B, 
Weiss et al. 31 who used a specially designed constant temperature anemometer to map out the freestream 
disturbance field within a Mach 2.54 short duration supersonic facility at Stuttgart University, and Masutti 
et al. 32 who characterized the freestream fluctuations in the Mach 6 Hypersonic Wind Tunnel H3 of the von 
Karrnan Institute. 

Table 4 and Table 5 list the freestream values of several fluctuating flow variables for Cases M6Tw025 and 
M6Tw076, respectively. Here, s is specific entropy, D is the vorticity, and the subscript ‘t’ denotes stagnation 
quantities. A comparison of the data from these two tables indicates that the magnitude of freestream 
fluctuations normalized by the respective mean values increase significantly as the wall temperature decreases, 
including both the velocity fluctuations and the fluctuations in thermodynamic variables. In particular, the 
pressure fluctuations in the free stream, including p' rms /p and p't rm s/Pti are considerably different for both 
cases (3.52% vs 2.01% for p' rms /p , 1.13% vs 0.648% for p' trms /p t , larger value for the cold-wall case). 
However, p' rms /p and p'trms/Pt bear nearly the same ratio of about 1.8 across the two cases. For both wall- 
temperature cases, the variation in rms amplitudes of velocity fluctuations along the three Cartesian axes 
indicates the anisotropy of the freestream velocity fluctuations, with the wall-normal velocity fluctuations 
the largest in the freestream. The relative perturbations in thermodynamic variables are nearly an order of 
magnitude larger than the velocity fluctuations. Consequently, the mass-flux fluctuations in the freestream 
are dominated by the density fluctuations, similar to the findings in the Mach 2.5 supersonic boundary 
layer. 14 Also observe that u! and p ' are not perfectly anti-correlated for either case, which is assumed by 
Laufer 5 in analyzing his hot-wire data. The negative correlation between u' and p' becomes stronger in the 
cold-wall case (—0.845 vs —0.632). 

The fact that s' rms /R « p' rms /p confirms the negligible contribution of entropy mode to the overall 
freestream disturbance. The pressure fluctuation in the freestream is nearly perfectly correlated with the 
fluctuations in density and temperature, so that the fluctuations in thermodynamic quantities satisfy the 
following isentropic relations 

V 

rms 

T 

Prms 

P 

analogous to purely acoustic disturbances. The correlation between the fluctuations in pressure and velocity is 
also strong, but significantly less than unity. The near zero correlation between pressure and spanwise velocity 
fluctuations is a consequence of spanwise homogeneity and simply confirms the satisfactory convergence of 
the statistical estimates. The relative importance of the acoustic and vorticity modes in the freestream 
can be assessed via the ratio of the dilatation variance ( dui/dx.i ) /2 and the vorticity variance fl'fl' in the 
freestream. These variances are representative of the disturbance energy in the acoustic and vorticity modes, 
respectively, and the large value of (dui/dxi)' 2 /Q,' i Ct' i implies the overwhelming dominance of acoustic mode 
over the vorticity mode in the free stream. 

The dominance of the acoustic mode over the other two modes confirms that a purely acoustic field in the 
free stream is successfully isolated by the present DNS. In typical high-speed (noisy) wind tunnels, however, 
freestream disturbances may originate from multiple disturbances sources (screens, heat exchangers, valves, 
bends, diffusers, and compressors upstream of the test section and vibration, roughness and waviness, bound- 
ary layer transition, and model- wall interference within the test section) and while acoustic disturbances may 
still be dominant, all three modes contribute to the freestream fluctuations in general. 2, 31-33 In that regard, 
the present simulation provides a more controlled setting for studying the spectral features and production 
mechanisms specific to the acoustic disturbance environment resulting from turbulent tunnel wall boundary 
layers. 

C. Frequency spectra of pressure fluctuations 

Figure 9 shows the comparison of frequency spectra at the wall and in the free stream for Cases M6Tw025 and 
M6Tw076. The wall-pressure spectrum shows a strong wall temperature dependence, which scales neither 
in outer variables (Figure 9a) nor in inner variables (Figure 9a). Figure 10a shows that the pre-multiplied 
wall-pressure spectrum for Case M6Tw025 consists of a sharper peak with a faster roll-off at high frequencies 
compared with Case M6Tw076, and wall cooling causes a shift in the dominant frequency from uj6/Uoo ss 9 
for Case M6Tw076 to uiS/U^ ~ 15 for Case M6Tw025. 
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For the freestream pressure spectrum, Figure 9c shows that the low-frequency range of <j> p (u ;) scales in 
outer variables; and Figure 9d shows that the high-frequency range scales in inner-variable scales, which 
conforms to the findings of wall-pressure spectrum in low-speed adiabatic flows. 34 Moreover, Figure 10b 
shows that the dominant frequency of freestream pressure spectrum is independent of wall temperature, 
indicating insignificant influence of wall cooling on the freestream pressure spectrum. 


D. Space-time correlation and propagation speed 


The statistical properties of the pressure field are investigated through the space-time correlation coefficient 
defined as 


Cpp{ Ax, Ay, At) 


p'(x, y, z, t)p'{x + Ax, y + Ay, z,t + At) 
(p' 2 (x,y,z,t)j 7 (p ,2 (x + Ax,y + Ay,z,t + At)) 


1/2 


( 4 ) 


where Aa; and Ay are spatial separations in the streamwise and spanwise directions, respectively, and At is 
the time delay. 

The space-time correlation contours C pp (Ax, 0, At) of the surface and freestream pressure fluctuations 
are shown in figures 11a and lib, respectively, for Cases M6Tw025 and M6Tw076. The skewed shape of 
the contours at both locations indicates the convective nature of the pressure held, which is characterized 
by downstream propagation of either the coherent pressure-carrying eddies within the boundary layer or the 
pressure wavepackets in the free stream. The similar inclination of the space-time correlation contours in 
the free stream between Cases M6Tw025 and M6Tw076 indicates that wall cooling has almost no influence 
on the overall propagation speed of radiated pressure wavepackets in the free stream. 

To quantitatively represent the overall propagation speed of pressure-carrying eddies or wavepackets, 
the bulk propagating speed (Ub) is defined as the one that minimizes the difference between the real time 
evolution of p(x,t) and a frozen wave p(x — Ubt). Following this definition, the following expression can be 

obtained 

TT = (dp/dt)(dp/dx) ^ 

b {. dp/dx ) 2 


The same definition of bulk propagation speed was introduced by Del Alamo 35 for the streamwise velocity 
fluctuations. 

Figure 12 plots the bulk propagation speed of the pressure fluctuation as a function of wall-normal 
distance. It is shown that wall cooling decreases the convection speed of pressure-carrying eddies within the 
boundary layer but has nearly no influence on the propagation speed of radiated pressure wavepackets in 
the free stream, where the bulk propagation speed is approximately OTC/qo for both wall-temperature cases. 

Figure 12 further shows that for both cases, the bulk propagation speed of the pressure fluctuation is 
larger than the local mean velocity near the wall and becomes smaller than the mean velocity farther away 
from the wall, which is consistent with previous findings. 13 ’ 14,27,36 The convection speed of pressure-carrying 
eddies is larger than the local mean velocity over a significant portion of boundary layer (up to z/S ss 0.65) 
for Case M6Tw076 while the analogous region is much smaller for the colder wall case (Case M6Tw025). 
The region where u ss Ub seems to be limited to the upper buffer layer and the log layer, at least for the 
present study. 

The value of convection/propagation speed is known to vary depending on specific definitions. 3 ' Thus, 
further comparison with different definitions of the convection speed over a wider range of flow conditions 
is required to confirm the above observations. Also, it is well known that the propagation/convection speed 
of turbulent eddies is a function of temporal frequency and spatial wave numbers; 35, 36 and further study is 
necessary to characterize the influence of wall cooling on the frequency- and wave-number dependence of the 
convection speed. 

Figures 13a and 13b show the variation of bulk convection speed with freestream Mach number at the 
wall and in the free stream, respectively. While the data by Kistler & Chen 38 shows that the convection 
speed at the wall plateaus for > 2, the current DNS and the data reported by Bernadini and Pirozzoli 27 
and Duan et al. 14 show a weak increase with the freestream Mach number. Both the current DNS and the 
experimental measurements by Laufer 5 show that the bulk convection speed in the free stream increases with 
freestream Mach number. All the freestream bulk convection speeds fall within the region where M r > 1, 
with M r = ([Zoo — Ub)/aoo, which supports the basic concept of ‘eddy Mach wave’ radiation 10, 11 and shows 
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that the Mach-wave-type radiation is produced by eddies, which are convected supersonically with respect 
to the free stream. 


E. Wave-front orientation 


In this section, the wave- front orientation of Mach-wave radiation is discussed, which is an important quantity 
from the standpoint of transition analysis, since the receptivity characteristics are known to be sensitive to 
the orientation of the plane-wave disturbance. 39 

Figure 14 shows that the instantaneous pressure field in the free stream for Cases M6Tw025 and 
M6Tw076. For both wall temperatures, the acoustic radiation consists of randomly spaced wavefronts, each 
with a limited spatial coherence. The wave fronts exhibit a preferred orientation within the streamwise- 
wall-normal (x-z) plane. By assuming purely plane-wave like radiation pattern, the wave-front orientation 
can be derived using the following acoustic relation 40 


< = 1 p' 

Uqo iMoo Poo 


( 6 ) 


where u n = u ■ n is the velocity normal to the wave front and n = (cos9 n , sin9 n ) is the plane- wave normal 
direction. The wave-front orientation that is statistically most likely can be determined to be the direction 
that minimizes the difference between u' n rma /Uoo and p' rms /{'yM 00 p 00 ). By using the freestream statistics 
analogous to Tables 5 and 4), 9 n w 120° (correspondingly the angle between the wave front and the flow 
direction 9 = 30°) is obtained for both cases, which agrees with the approximate inclination of the wave 
fronts from the numerical schlieren image in figure 14. Laufer 5 has used a similar but less rigorous relation 
v! n cos9 n = u' to estimate the wave-orientation, since the streamwise velocity fluctuation u' is the only velocity 
component that could be measured in his experiments. 

The reader should note that the freestream acoustic field does not correspond to truly planar waves as 
indicated by the instantaneous pressure field in the spanwise-wall-normal (y-z) plane for both cases. The 
deviation from purely planar behavior is also evident in a statistical sense from the imperfect correlation 
between p' and the streamwise ( v ! ) and wall-normal ( w ') velocity fluctuations in the freestream region. 
(Recall the data presented in Tables 5 and 4). The finite spanwise extent of the pressure wavepackets is 
consistent with the finite size of acoustic sources that are responsible for generating the waves. 

While the overall orientation of the acoustic radiation for Cases M6Tw025 and M6Tw076 is nearly the 
same (9 ss 30°), the wavepackets of the radiation field for Case M6Tw025 are more spatially coherent 
compared with Case M6Tw076, as indicated by the larger spatial extent of the radiated wavefronts in the 
numerical schlieren image. Such an increase in spatial coherence for Case M6Tw025 is consistent with the 
larger correlations between p' and (u')~ and (u/)-velocity fluctuations. The increased coherence of radiated 
wavepackets with wall cooling may be a result of increased coherence of near-wall turbulent structures due to 
wall cooling, 17 as the acoustic sources that can directly radiate to the free stream are associated with eddies 
with M r > 1, and, hence, are confined to the near- wall region of the boundary layer. Detailed analysis of 
the structural characteristics of the acoustic sources is still necessary to confirm this argument. 


IV. Conclusion 

DNS of Mach 6 turbulent boundary layers with two different wall temperatures ( T w /T r = 0.25, 0.76) are 
conducted to investigate the effect of wall temperature on the pressure fluctuations generated by hypersonic 
boundary layers. Simulations show that wall cooling significantly modifies the pressure-fluctuation intensities 
near the wall, with p' w ^ rms /T w varying from 2.6 for T w /T r = 0.76 to 3.4 for T w /T r = 0.25. Furthermore, the 
frequency spectra of wall-pressure fluctuations for the two cases show considerable differences when plotted 
in terms of either outer-layer or inner-layer variables. The peak of the pre-multiplied spectrum shifts to a 
higher value as the wall temperature decreases. 

In terms of the freestream pressure fluctuations, although the radiation intensity shows a strong wall- 
temperature dependence when normalized by the mean freestream pressure (p^), the radiation intensity 
compares well between the two cases when normalized by the local wall shear. Wall cooling increases the 
spatial correlation of the freestream pressure field, but has little influence on the frequency spectrum, the 
propagation speed, or the dominant orientation of the radiation wave fronts. 


8 of 21 


American Institute of Aeronautics and Astronautics 



Acknowledgments 


This material is based on the work supported by National Aeronautics and Space Administration, Langley 
Research Center under the Research Cooperative Agreement No. NNL09AA00A (subcontracted through the 
National Institute of Aerospace). The authors would like to thank Dr. Katya Casper for providing her wind- 
tunnel measurements for comparison with DNS. Prof. Pino Martin of the University of Maryland provided the 
original code which has been modified for the current study. The simulations have been conducted using the 
Pleiades Supercomputer of the NASA Advanced Supercomputing (NAS) Division and the Extreme Science 
and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant 
number OCI-1053575. 


References 

1 Pate, S. R. and Schuller, C. J., “Radiated Aerodynamic Noise Effects on Boundary Layer Transition in Supersonic and 
Hypersonic Wind Tunnels,” AIAA Journal , Vol. 7, No. 3, 1969, pp. 450-457. 

2 Donaldson, J. and Coulter, S., “A Review of Free- Stream Flow Fluctuation and Steady-State Flow Quality Measurements 
in the AEDC/VKF Supersonic Tunnel A and Hypersonic Tunnel B,” AIAA Paper 95-6137, 1995. 

3 Schneider, S. P., “Effects of High-Speed Tunnel Noise on Laminar-Turbulent Transition,” Journal of Spacecraft and 
Rockets , Vol. 38, No. 3, 2001, pp. 323—333. 

4 Borg, M. P. and Schneider, S. P., “Effect of Free-stream Noise on Roughness-Induced Transition for the X-51A Forebody,” 
Journal of Spacecraft and Rockets , Vol. 45, No. 6, 2008, pp. 1106-1116. 

5 Laufer, J., “Some Statistical Properties of the Pressure Field Radiated by a Turbulent Boundary Layer,” Physics of 
Fluids, Vol. 7, No. 8, 1964, pp. 1191-1197. 

6 Stetson, K. F., “Nosetip Bluntness Effects on Cone Frustrum Boundary-Layer Transition in Hypersonic Flow,” AIAA 
Paper 83-1763, 1983. 

7 Bushnell, D. M., “Notes on Initial Disturbance Fields for the Transition Problem,” Instability and Transition, edited by 
M. Y. Hussaini and R. G. Voigt, Springer- Verlag, Berlin, Vol. 1, 1990, pp. 217-232. 

8 Kendall, J. M., “Supersonic Boundary Layer Transition Studies,” Space Program Summary , Vol. 3, 1970, pp. 43-47. 

9 Beresh, S. J., Henfling, J. F., Spillers, R. W., and Pruett, B. O. M., “Fluctuating Wall Pressures Measured beneath a 
Supersonic Turbulent Boundary Layer,” Physics of Fluids, Vol. 23, No. 7, 2011, 075110. 

10 Phillips, O. M., “On the Generation of Sound by Supersonic Turbulent Shear Layers,” Journal of Fluid Mechanics, 
Vol. 9, 1960, pp. 1-28. 

11 Ffowcs- Williams, J. E. and Maidanik, G., “The Mach Wave Field Radiated by Supersonic Turbulent Shear Flows,” 
Journal of Fluid Mechanics, Vol. 21, 1965, pp. 641-657. 

12 Duan, L., Choudhari, M., and Wu, M., “Numerical Study of Pressure Fluctuations Due to High-Speed Turbulent Bound- 
ary Layers,” AIAA Paper 2012-3070, 2012. 

13 Duan, L. and Choudhari, M. M., “Numerical Study of Pressure Fluctuations due to a Mach 6 Turbulent Boundary 
Layer,” AIAA Paper 2013-0532, 2013. 

14 Duan, L., Choudhari, M. M., and Wu, M., “Numerical Study of Acoustic Radiation due to a Supersonic Turbulent 
Boundary Layer,” Journal of Fluid Mechanics, Vol. 746, 2014, pp. 165-192. 

15 Welch, P. D., “The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time 
Averaging Over Short, Modified Periodograms,” IEEE Trans. Audio Electroacoustics, Vol. AU-15, 1967, pp. 70-73. 

16 Martfn, M., “DNS of Hypersonic Turbulent Boundary Layers. Part I: Initialization and Comparison with Experiments,” 
Journal of Fluid Mechanics, Vol. 570, 2007, pp. 347-364. 

17 Duan, L., Beekman, I., and Martin, M. P., “Direct Numerical Simulation of Hypersonic Turbulent Boundary Layers. Part 
2: Effect of Wall Temperature,” Journal of Fluid Mechanics , Vol. 655, 2010, pp. 419-445. 

18 Duan, L., Beekman, I., and Martin, M. P., “Direct Numerical Simulation of Hypersonic Turbulent Boundary Layers. Part 
3: Effect of Mach Number,” Journal of Fluid Mechanics, Vol. 672, 2011, pp. 245-267. 

19 Duan, L. and Martin, M. P., “Direct Numerical Simulation of Hypersonic Turbulent Boundary Layers. Part 4: Effect of 
High Enthalpy,” Journal of Fluid Mechanics, Vol. 684, 2011, pp. 25—59. 

20 Wu, M. and Martin, M. P., “Direct Numerical Simulation of Supersonic Boundary Layer over a Compression Ramp,” 
AIAA Journal, Vol. 45, No. 4, 2007, pp. 879-889. 

21 Wu, M. and Martin, M. P., “Analysis of Shock Motion in Shock Wave and Turbulent Boundary Layer Interaction using 
Direct Numerical Simulation Data,” Journal of Fluid Mechanics, Vol. 594, 2008, pp. 71-83. 

22 Priebe, S. and Martin, M. P., “Low-Frequency Unsteadiness in Shock Wave Turbulent Boundary Layer Interaction,” 
Journal of Fluid Mechanics, Vol. 699, 2012, pp. 1-49. 

23 Taylor, E. M., Wu, M., and Martin, M. P., “Optimization of Nonlinear Error Sources for Weighted Non- Oscillatory 
Methods in Direct Numerical Simulations of Compressible Turbulence,” Journal of Computational Physics, Vol. 223, No. 1, 
2006, pp. 384-397. 

24 Harris, J. and Blanchard, D., “Computer Program for Solving Laminar, Transitional, or Turbulent Compressible 
Boundary-Layer Equations for Two-Dimensional and Axisymmetric Flow,” NASA-TM-83207, 1982. 

25 Bookey, P., Wyckham, C., Smits, A. J., and Martin, M. P., “New Experimental Data of STBLI at DNS/LES Accessible 
Reynolds Numbers,” AIAA Paper 2005-309, 2005. 


9 of 21 


American Institute of Aeronautics and Astronautics 



26 Pirozzoli, S. and Bernardini, M., “Turbulence in Supersonic Boundary Layers at Moderate Reynolds Numbers,” Journal 
of Fluid Mechanics , Vol. 688, 2011, pp. 120-168. 

27 Bernardini, M. and Pirozzoli, S., “Wall Pressure Fluctuations beneath Supersonic Turbulent Boundary Layers,” Physics 
of Fluids, Vol. 23, No. 8, 2011, 085102. 

28 Kovasznay, L. S. G., “Turbulence in Supersonic Flow,” Journal of Aeronautical Sciences , Vol. 20, 1953, pp. 657-674. 

29 Logan, P., “Modal Analysis of Hot-Wire Measurements in Supersonic Turbulence,” AIAA Paper 88-423, 1988. 

30 Smits, A. J. and Dussauge, J. P., Turbulent Shear Layers in Supersonic Flow, American Institute of Physics, 2nd ed., 
2006. 

31 Weiss, J., Knauss, H., and Wagner, S., “Experimental Determination of the Freestream Disturbance Field in a Short- 
Duration Supersonic Wind Tunnel,” Experiments in Fluids, Vol. 35, 2003, pp. 291-302. 

32 Masutti, D., Spinosa, E., Chazot, O., and Carbonaro, M., “Disturbance Level Characterization of a Hypersonic Blow- 
Down Facility,” AIAA Journal, Vol. 50, No. 12, 2013, pp. 2720-2730. 

33 Morkovin, M. V., “On Transition Experiments at Moderate Supersonic Speeds,” Journal of the Aeronautical Sciences, 
Vol. 24, No. 7, 1957, pp. 480-486. 

34 Bull, M. K., “Wall-Pressure Fluctuations beneath Turbulent Boundary Layers: Some Reflection on Forty Years of Re- 
search,” Journal of Sound and Vibration, Vol. 190, No. 3, 1996, pp. 299-315. 

35 Del Alamo, J. C. and Jimenez, J., “Estimation of Turbulent Convection Velocities and Corrections to Taylor’s Approxi- 
mation,” Journal of Fluid Mechanics, Vol. 640, 2009, pp. 5-26. 

36 Kim, J. and Hussain, F., “Propagation Velocity of Perturbations in Turbulent Channel Flow,” Physics of Fluids, Vol. 5, 
No. 3, 1993, pp. 695-706. 

37 Choi, H. and Moin, P., “On the Space-Time Characteristics of Wall-Pressure Fluctuations,” Physics of Fluids, Vol. 2, 
No. 8, 1990, pp. 1450-1460. 

38 Kistler, A. L. and Chen, W. S., “The Fluctuating Pressure Field in a Supersonic Turbulent Boundary Layer,” Journal 
of Fluid Mechanics, Vol. 16, 1963, pp. 41-64. 

39 Fedorov, A. V., “Receptivity of a High-Speed Boundary Layer to Acoustic Disturbances,” Journal of Fluid Mechanics, 
Vol. 491, 2003, pp. 101-129. 

40 Liepmann, H. W. and Roshko, A., Elements of Gasdynamics, John Wiley Sz Sons, Inc. New York, 1957. 


10 of 21 


American Institute of Aeronautics and Astronautics 



Table 1. Freestream conditions for Mach 6 DNS of turbulent boundary layers. 


/Too (m/s) Poo(kg/m 3 ) T^K) 

5.86 870.4 0.0427 55.0 


Table 2. Boundary layer properties at the station selected for the analysis of the pressure field 
( x a = 57(5, for Case M6Tw076 and x a = 80(5; for Case M6Tw025). The local boundary layer 
thickness 8 is approximately 5 « 1.8(5; for Case M6Tw076 and 8 « 2.6(5; for Case M6Tw025. 


Case 

^iu(K) 

T w /T r 

Re e 

Re T 

Res2 

6 (mm) 

H 

S (mm) 

z T (fim) 

u T (m/s) 

M6Tw025 

97.5 

0.25 

4650.7 

465.8 

2347.1 

0.193 

8.45 

3.70 

7.98 

34.0 

M6Tw076 

300 

0.76 

9656.4 

464.1 

1782.8 

0.969 

13.6 

24.4 

52.6 

45.0 


Table 3. Grid resolution and domain size for Case M6Tw076. L x , L y , and L z are the domain 
size in the streamwise, spanwise, and wall-normal directions, respectively. Ax + and A y + are the 
uniform grid spacing in the streamwise and spanwise directions, respectively, with the viscous 
length scale z T = 52.6/im corresponding to the station selected for statistical analysis (x/Si = 57). 
A 2 + irl and A 4m are the minimum and maximum wall-normal grid spacing for 0 < z/8i < 7. 


Case 

N x x N y x N z 

Lx/5, 

Ly/Si 

L z /5i 

Ax + 

Ay+ 

Az Ln 

A^+ 

“ max 

M6Tw076, Baseline 

1600 x 800 x 500 

58.7 

15.7 

39.7 

9.64 

5.14 

0.51 

5.33 

M6Tw076, AI 

1920 x 320 x 500 

70.4 

6.26 

39.7 

9.64 

5.14 

0.51 

5.33 

M6Tw076, All 

2400 x 480 x 700 

58.7 

6.26 

39.7 

6.43 

3.43 

0.51 

3.55 
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Table 4. The disturbance field in the free stream for Case M6Tw025. 


^ rms ^ rms / ^ Prms/P Prms/ P r ^rms/ r ^' 

2.45 x 10 -3 1.64 x 10" 3 3.26 x 1CT 3 3.53 x 1CT 2 2.52 x 10" 2 1.01 x 10" 2 

(ouY hm T! IT * J /o, OuJdxA ' 2 Ifl'Sl'- s' JR 


(, P U )'rms/P U T 't ,rms /^t Pt,rms/Pt 

2.31 x 1CT 2 3.23 x 1CT 3 1.13 x 10" 2 


(du i /dx i y 2 /n , i p.' i s' rm JR 
2985 2.18 x 1CT 4 


IL p / U rrns p rrns V p /‘U rrns P rrns W P / ^XrmsPrms P P / PrmsPrms ^ P /^rmsPri 


-0.845 


-0.039 


0.955 


1.0 


1.0 


Table 5. The disturbance field in the free stream for Case M6Tw076. 


'"’rms/'"’ ^ rms /'"' ^ rms / ^ Prms/P Prms/ P T'rms /'^ 1 

1.31 x 1(T 3 1.04 x 10" 3 2.05 x 1CT 3 2.01 x 1CT 2 1.43 x 10" 2 5.78 x 1CT 3 

(/ PU)'rms/W Tirms/Tt P't,rms/Pt {dUi/dXi)' 2 /fi'ft' s' rms /R 

1.35 x KT 2 1.93 x 10" 3 6.48 x 10" 3 1131.2 2.27 x 1CT 3 


U p /u rrns P rms X P / X rrns P rrns W p / W rrns P rrns P P /PrmsPrms 


T' p' /T' rrns p' rl 


— 0.632 


— 0.005 


0.918 


0.999 


0.993 



Figure 1. Computational domain and simulation setup for Case M6Tw025, which is also repre- 
sentative of the other case. The reference length <5; is the thickness of the boundary layer (based 
on 99% of the freestream velocity) at the inflow plane. An instantaneous flow field is shown in 
the domain, visualized by an iso-surface of the magnitude of density gradient, |Vp|<5i/pc>o = 0.98, 
colored by the streamwise velocity component (with levels from 0 to [Zoo, blue to red). 
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Figure 2. The correlation coefficient of pressure fluctuations as a function of spanwise separation 
at the wall (z/5 = 0) and in the freestream (z/<5 = 2.43), at selected streamwise locations x = 60 Si 
and x = 85<5i for Case M6Tw025. 




Figure 3. (a) Streamwise and (6) wall-normal distributions of the normalized rms pressure 

fluctuation in the downstream portion of the computational domain (60 < x/8i < 85) for Case 
M6Tw025. The freestream location is picked as z/8 = 2.43. 
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Figure 4. Comparison of (a) rms pressure fluctuation and ( b,c ) frequency spectra of pressure 
signal at the wall and in the freestream for DNS solutions based on varying domain size and grid 
resolutions. The DNS cases are listed in Table 3. 
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Figure 5. Comparison of DNS results (Case M6Tw076) with those of a Mach-5.8 turbulent 
boundary layer on the nozzle wall of the Boeing/AFOSR Mach-6 Quiet Tunnel under noisy-flow 
conditions {Re = 9.69 x 10 6 /m, Pt,oo = 965 kPa, Tt,oo = 429 K). The wind-tunnel measurement and 
the calculation using Flarris boundary- layer Code 24 were conducted by Katya Casper, (o) Mean 
velocity profile; ( b ) Mach number profile; (c) frequency spectrum in outer scale; ( d , ) frequency 
spectrum in inner scale. 
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Figure 6. (a) van Driest transformed mean velocity profile (k = 0 . 41 , C = 5 . 5 ) and ( b ) van Driest 
transformed mean deficit velocity for Case M6Tw025. Symbols denote the DNS by Pirozzoli & 
Bernardini 26 at Mach 2, Re T = 497 (squares) and the experiment by Bookey et al. 25 at Mach 2.9, 
Re T = 501 (down-pointing triangles). The DNS cases are listed in Table 2. 
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Figure 7. (a) Pressure fluctuation rms normalized by the wall-shear stress for the present DNS 
and the Mach 2.5 DNS by Duan et al. 14 in outer and inner scales. 
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Figure 8. Intensity of freestream pressure fluctuation compared with the experiments by Laufer. 5 
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Figure 9. Comparison of power spectral density at the wall and in the freestream (z/5 = 2.43) for 
Cases M6Tw025 and M6Tw076. 
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Figure 10. Comparison of pre-multiplied power spectral density at the wall and in the freestream 
(z/5 = 2.43) for Cases M6Tw025 and M6Tw076. 



Figure 11. (a) space-time correlation at the wall and in the freestream (z/5 = 2.43). The solid 
contours are for Case M6Tw025 and the dashed contours are for Case M6Tw076. 
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Figure 12. Bulk propagation speed of pressure fluctuations as a function of wall-normal distance 
in (a) outer and (b) inner units. 
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Figure 13. Bulk convection speeds of the pressure fluctuation as a function of freestream Mach 
number: (a) at the wall; (b) in the freestream. Symbols: squares, Kistler & Chen; 38 Left triangles, 
Bernardini & Pirozzoli; 27 diamonds, Laufer; 5 deltas, Duan et al.; 14 right triangles, Present DNS, 
Case M6Tw025; gradients, Present DNS, Case M6Tw076. Lines: M r = 1. 
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(b) M6Tw076 


Figure 14. Numerical schlieren image based on instantaneous flow field for Cases M6Tw025 and 
M6Tw076. Contour levels are selected to emphasize disturbances in the freestream. 6 is the 
angle between the wave front and the flow direction; and the vertical dashed line indicates the 
streamwise location of the selected spanwise— wall-normal plane visualized in the right panel. 
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